In [7]:
emissionssw <- read.table("datasets/emissionssw.dat",header=TRUE)
n <- nrow(emissionssw)
In [8]:
summary(emissionssw)
nox noxem ws humidity Min. : 1.2 Min. : 102.5 Min. : 0.100 Min. :0.009332 1st Qu.: 47.7 1st Qu.: 675.0 1st Qu.: 1.033 1st Qu.:0.060673 Median : 88.3 Median :2183.6 Median : 1.667 Median :0.100842 Mean :116.0 Mean :2250.7 Mean : 2.087 Mean :0.171872 3rd Qu.:152.5 3rd Qu.:3741.7 3rd Qu.: 2.750 3rd Qu.:0.168508 Max. :694.5 Max. :5362.1 Max. :11.367 Max. :4.099019
In [118]:
with(options(repr.plot.width = 12, repr.plot.height = 8), {
par(pin = c(20, 20), mfrow = c(4, 1), mar = c(2.5, 2, 1.5, 1), mgp = c(1.5, 0.7, 0))
par(cex.lab = 1, cex.axis = 1, cex.main = 1.5, cex.sub = 1.2)
plot(emissionssw$nox, type = 'l', main = 'NOx (PPB)', xlab = '', ylab = '')
plot(emissionssw$noxem, type = 'l', main = 'NOx Emissions (unknown)', xlab = '', ylab = '')
plot(emissionssw$ws, type = 'l', main = 'Wind Speed (m/s)', xlab = '', ylab = '')
plot(emissionssw$humidity, type = 'l', main = 'Humidity (g/kg)', xlab = 'Index', ylab = '')
})
In [9]:
cor(emissionssw)
| nox | noxem | ws | humidity | |
|---|---|---|---|---|
| nox | 1.00000000 | 0.51758582 | -0.34754186 | -0.04739774 |
| noxem | 0.51758582 | 1.00000000 | 0.11653505 | -0.03381921 |
| ws | -0.34754186 | 0.11653505 | 1.00000000 | 0.01722212 |
| humidity | -0.04739774 | -0.03381921 | 0.01722212 | 1.00000000 |
In [98]:
with(options(repr.plot.width = 8, repr.plot.height = 8), {pairs(emissionssw)})